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ABSTRACT 

We present a Faraday rotation measure (RM) study of the diffuse, polarized, radio emission from the 
giant lobes of the nearest radio galaxy, Centaurus A. After removal of the smooth Galactic foreground 
RM component, using an ensemble of background source RMs located outside the giant lobes, we are 
left with a residual RM signal associated with the giant lobes. We find the most likely origin of this 
residual RM is from thermal material mixed throughout the relativistic lobe plasma. The alternative 
possibility of a thin-skin/boundary layer of magnetoionic material swept up by the expansion of the 
lobes is highly unlikely since it requires, at least, an order of magnitude enhancement of the swept up 
gas over the expected intragroup density on these scales. Strong depolarisation observed from 2.3 to 
0.96 GHz also supports the presence of a significant amount of thermal gas within the lobes; although 
depolarisation solely due to RM fluctuations in a foreground Faraday screen on scales smaller than 
the beam cannot be ruled out. Considering the internal Faraday rotation scenario, we find a thermal 
gas number density of ~ ICE 4 cm~ 3 implying a total gas mass of ~ 10 10 M Q within the lobes. The 
thermal pressure associated with this gas (with temperature kT ~ 0.5 keV, obtained from recent 
X-ray results) is approximately equal to the non-thermal pressure, indicating that over the volume of 
the lobes, there is approximate equipartition between the thermal gas, radio-emitting electrons and 
magnetic field (and potentially any relativistic protons present). 

Subject headings: radio galaxies: individual: Centaurus A (NGC 5128) — radio galaxies: magnetic 
fields 



1. INTRODUCTION 

The integrated history of a radio galaxy is encoded in 
the magnetized, relativistic gas within its lobes. These 
lobes have been inflated by outflows emanating from 
the cent ral supermassive black hole (SMBH) of the host 
galaxy (Begelman et al. 1984) with their structure illu- 
minated by the synchrotron emission we detect in the 
radio band. Through their formation, radio galaxies can 
deposit 10 55 -10 62 erg of mechanical energy into the in- 
tergalactic medium (IGM ) over the lifetime of the source 
( |McNamara et aL||2009 ) , as we ll as potentially enriching 
the IGM with heavy elements ( |Aguirre jrtjdj[2001 |Rj3U- 
l land et al.|2007| ) and magnetic helds ( |Furlanetto & Loeb 
|2U0l| . Determining exactly how radio galaxies transfer 
energy and material to the IGM is essential for under- 
standing their impact on both ga l axy and cosmic struc- 
ture evolution jBower et al ||2006| |Best et al.||2006| |Cro- 
ton et alT||2006 1 . In particular, knowledge ot the traction 
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Faraday rotation provides a sensitive diagnostic for the 



ot "non" -radiating particles (thermal gas and relativis 
tic protons) is of primary importance for understanding 
the pressure balance and dynamics of radio galaxy lobes 



with respect to their environments jMorganti et alTl988 
Dunn et al.|2005l|Birzan et al.|2008"f|Croston et al |2008 ) 



as well as their potential to provide seed particles to r high 
energy emission and ultra-high energy cosmic rays ( Abdo 
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et al. 2010). 
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thermal gas within radio galaxy lobes of < 10~ 4 cm~ 3 , 
assuming the minimum-pressure magnetic field intensity. 
However, conclusively disentangling the different contri- 
butions from Faraday rotation that is internal and/or 
external to radio galaxy lobes has proven very difficult. 

In this paper, we present complementary spectropolari- 
metric observations of the nearest radio galaxy Centau- 
rus A at 1.4 GHz with both the Parkes 64 m single-dish 
radio telescope and the Australia Telescope Compact Ar- 
ray (ATCA) aperture synthesis telescope. The Parkes 
observations are sensitive to diffuse emission which al- 
lows us to measure the Faraday rotation measure (RM) 
of the polarized emission from Centaurus A itself, while 
the h igh an gular resolution ATCA observations ( |Feain| 



s.o'sulli van@physics.usyd.edu.au. 



et al. 2009) provide us with an ensemble of RMs from 
281 background radio sources that are uncontaminated 
by any diffuse polarized emission and are located along 
sightlines both inside and outside the giant lobes. Cen- 
taurus A is the only radio galaxy for which such an anal- 
ysis is possible with current sensitivity due to its large 
angular size of ~9° x 4°. For a di stance to Centaurus A 
of ^3.8 Mpc (Harris et al. 2010), 1° corresponds to a 
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projected linear size of ~66 kpc. 

In Section 2, we describe the observations, calibration 
and RM analysis. Section 3 outlines our results along 
with a discussion on various models and the implication 
for the giant lobes. Our conclusions are presented in 
Section 4. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Parkes Data 

A 10° x 14° area centred on the position of Centau- 
rus A was observed using the Parkes 64 m telescope at 
1.4 GHz over a period of 10 days (2009 March 2-11), 
to talling 80 hrs of te lescope time. The Parkes data used 



Feain et al. (2011 ) were archival but the new observa- 
tions presented here were required because the archival 
data did not contain polarisation information. The ob- 
servations were conducted using a cross-scanning tech- 
nique wit h the H-OH receiv er and a quarter-wave plate, 
similar to Mao et al. (2012), providing 0.25 MHz chan- 
nels across 256 MHz of bandwidth centerd on 1388 MHz. 
The total useable bandwidth covered 1312 to 1480 MHz. 
The FWHM for the Parkes telescope at 1380 MHz is ap- 
proximately 14.4'. Data were obtained while scanning 
the telescope at a rate of 3°. 5 min -1 and sampling ev- 
ery 1 s, meaning that data were recorded every 7'. This 
resulted in 120 right ascension (RA) scans and 86 dec- 
lination (Dec) scans. The data were then imaged with 
8' x 8' pixels to ensure at least one data sample per pixel 
(some blank pixels occurred towards the Northern edge 
of the image where the scanning rate was slightly too 
fast for 7' x 7' pixels). This effectively smooths the im- 
age and we measure an actual gridded beam width of 
~ 16' at 1312 MHz. 

Calibration of all scans was done using the Parkes Con- 
tinuum Polarimetry Software (ParkesPol|^| package. 
The source PKS B1934 — 638 was used for flux and band- 
pass calibration due to its bright, s table and well k nown 
flux level of 14.95 Jy at 1380 MHz (|Reynolds|1994[ ). Ba- 
sic flagging was done on all scans before the 0.25 MHz 
channels were rebinned into 8 MHz channels and the po- 
larization calibration was performed on each channel sep- 
arately. Ten scans of 3C 138 at a range of feed angles 
from —45° to +45° provided a sufficient range of paral- 
lactic angle coverage to calibrate the instrumental polar- 
izations. The polarisation angle calibration was achieved 
using the known polarization angle of PKS B0043 — 424 
(+143.3° at 1.4 GHz, RM = +2 rad m -2 , Carretti & 
Haver korn, private communication). 

Stokes /, Q and U images were created for each fre- 
quency channel by combining all RA and Dec scans. Ini- 
tially, separate RA and Dec maps were created with a 
linear baseline removal performed using the edge of the 
image. This removes both the ground emission contam- 
ination as well as any structure on scales the size of the 
map and greater. The final maps are produced by com- 
bining the ortho gonal scan maps in Fourie r space using 



the technique of |Emerson fc Graeve (19881 which mini- 
mizes the effects of baseline drifts on raster scanned data. 

2.2. ATCA Data 
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Fig. 1. — Polarized intensity (in Jy beam -1 ) of the full observed 
region overlaid with the 1.4 GHz total intensity contours start- 
ing at 250 mjy beam -1 and increasing by factors of two. Diffuse 
polarized emission is detected throughout the entire region due to 
the Galactic foreground. The degree of polarization of the strongly 
polarized regions (p > 300 mjy beam -1 ) within the lobes ranges 
from ~ 10-40%. 

scope Compact Array (ATCA) at 1.4 GHz. They derived 
RMs for 281 background radio sources using essentially 
the same spectropolarimctric techniques as employed in 
this paper (see Section 2.3) with 24 x 8 MHz channels 
covering 1288 to 1480 MHz. In order to analyse the 
background radio sources without contamination from 
the large scale emission from Centaurus A, they filtered 
out the emission on large spatial scales by disregard- 
ing all visibilities from baselines shorter than 300 metres 
(1.4 kA). The high angular resolution ATCA observa- 
tions (~40") thus provided an ensemble of background 
radio sources that were uncontaminated by any diffuse 
polarized emission and located along sight-lines both in- 
side and outside the edge of the giant lobes. 

2.3. Rotation Measure Analysis 

Faraday rotation provides a direct diagnostic for the 
presence of magnetized thermal material though obser- 
vations of the change in the state of polarization with 
wavelength as the radiation passes through magnetoionic 
media on its path to us. The Faraday depth (0) of a par- 
ticular region of polarized emission is defined as 



<f>(L) = 0.81 / n e Budl rad m~ 



Feain et al. ( 2009 ) observed an ~45 deg 2 area centred 
on the host galaxy NGC 5128 with the Australia Tele- 

1 https:/ /svn. atnf.csiro.au/trac/parkespol 



(1) 



where n e is the electron number density in cm -3 , Bu is 
the line-of-sight magnetic field strength in jiG and L is 
the distance through the magnetoionic region in parsecs. 
The total observed Faraday rotation measure (RM), de- 
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Fig. 2. — Faraday rotation measure (RM) derived at each 8' X 8' 
pixel, overlaid with the 1.4 GHz total intensity contours. The range 
of RM on the color scale goes from —175 to +100 rad m -2 . 

fined as dx/dX 2 , where x is the polarisation angle, can 
have contributions from multiple regions along the line of 
sight. In the simplest case of a uniform external Faraday 
screen, the RM and Faraday depth are identical and the 
change in the polarization angle follows the relation 



X = Xo + <M , 



(2) 



where Xo is the intrinsic (A = 0) polarization angle. 

To determine the RM we first constructed cubes of 
Stokes Q and U vs. A 2 and applied the RM synthesis and 
RMC LEAN ( |Brentjens fc de Bruyn|[2005| fkeald et ah 
2009|) techniques follow ing the algorithm described by 



ullivan et al. (2012). This generated a cube of po- 



larized intensity (p = y/ Q 2 + U 2 ) at all Faraday depths 
from -3000 to +3000 rad m -2 from which the RM im- 
age was generated by extracting the RM at the peak 
polarized intensity. With the correct RM at each pixel 
the cube of polarized intensity is collapsed to form the 
polarized intensity image with the full sensitivity of the 
168 MHz band. 

Figure [T] shows the diffuse polarized emission detected 
from Centaurus A as well as from the Galaxy, which 
fills the entire field. The resolution in Faraday depth 
space (Scf) of our experiment, defined by the FWHM 
of the Rotation Measure Spread Function (RMSF), is 
310 rad uT 2 . The largest magnitude Faraday depth we 
can reliably detect is ^3000 rad m -2 and the maximum 
detectable Faraday thickness is ~80 rad m~ 2 . 

For a typical signal-to-noise ratio (S/N) of 50 to 100 
on source, the error in the RM estimate at each pixel 
is ~l-3 rad m" 2 (ARM ~ 8<j>/(2 x S/N)). This error 
estimate relies on the assumption that there is only one 
dominant Faraday depth component within the FWHM 
of the RMSF; if there are multiple Faraday depth compo- 



nents within our RMSF then we cannot uniquely specify 
their locations to better than 310 rad m~ 2 . However, for 
large S/N ratios the RMCLEAN algorithm may in prin- 
ciple be able to resolve components within the RMSF. 
We do not find any evidence for such components or for 
any significant broadening of the RMSF within the lobes 
of Centaurus A. In some regions of the observed field we 
find more than one peak in our Faraday depth spectrum. 
While the origin of these multiple Faraday depth compo- 
nents along the line of sight is interesting, none of these 
regions lie within the lobes of Centaurus A and are not 
included in our current study. 

3. RESULTS & DISCUSSION 

The measured RMs from the Parkes observations are 
shown in Figure [2] displaying a clear enhancement in the 
magnitude of the RM within the lobes of Centaurus A. 
Note that even with no Faraday rotating material in the 
vicinity of Centaurus A, we could observe such an en- 
hancement in the magnitude of the observed RM in re- 
gions where the polarized emission from Centaurus A 
dominates over the Galactic emission. This can occur be- 
cause the polarized emission of Centaurus A probes the 
full Faraday depth along the line of sight to us, while the 
diffuse Galactic polarized emission originates somewhere 
within our Galaxy and may only probe more local Fara- 
day rotating regions depending on the exact location(s) 
of the polarized emission within our Galaxy. However, 
due to potential changes in the line of sight magnetic 
field direction, a longer line of sight does not necessar- 
ily mean a larger observed Faraday depth. The Galactic 
emission may also be subject to strong internal Faraday 
rotation effects resulting in a significantly lower RM than 
that measured from the polarized emission from Centau- 
rus A. 

3.1. Residual RM from Centaurus A 

The mean RM of the background sources outside the 
lobes is —52.9 rad m -2 with a standard deviation of 
29.2 rad m" 2 (|Feain et al |2009| ). Therefore, the Faraday 
depth of our Galaxy is the dominant contributor to the 
total observed RM. The RMs from the 160 background 
sources located outside the lobes were used to remove the 
smooth part of the foreground Galactic Faraday rotating 
material. This was achieved by fitting a first-order two- 
dimensional polynomial to the point-source RMs outside 
the lobes. By subtracting this RM surface, with a gra- 
dient of ~ 6 rad m -2 deg -1 , from both the background 
source RMs (inside and outside the lobes) and the RM 
of Centaurus A, we obtain the residual RM signal. The 
colored pixels in Figure [3] show the residual RM signal 
from the polarized emission of Centaurus A while the 
open circles show the residual RM from the background 
sources inside the lobes. 

The mean of the magnitude of the residual RM from 
the polarized emission from Centaurus A greater than 
300 mJy beam' 1 is (|RM|) = 12.0 ± 0.3 rad nT 2 (us- 
ing the standard error of the mean^J The mean value 
of the RM from the Galactic polarized emission outside 
the lobes is approximately —27 rad m -2 with a standard 

2 Using only the same sightlines as the background sources we 
get <|RM|> = 13.0 ± 1.9 rad m" 2 . 
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deviation of 23 rad m -2 while the mean polarized inten- 
sity is ~ 52 mJy with a standard deviation of 20 mJy. 
As the Galactic polarized emission becomes a significant 
fraction of the Centaurus A polarized emission towards 
the edges of the lobes, the systematic error in the mea- 
sured RM will become large. Since we are unable to re- 
solve the two contributions with our current dataset, we 
have only included in our analysis the regions within the 
lobes where the polarized emission from Centaurus A is 
strong enough to make the effect from the Galactic emis- 
sion negligible (i.e., for Centaurus A polarized emission 
>300 mJy beam -1 , see Appendix). 

It is possible that the residual RM signal could be ex- 
plained by a degree-scale variation in the Galactic Fara- 
day depth that happens to coincide with the orientation 
of the giant lobes and is not sampled by the background 
sources outside the lobes. We consider this unlikely but 
cannot definitively rule out such a scenario. There is a 
clear asymmetry in the distribution of signs of the resid- 
ual RM between the Northern and Southern lobes. How- 
ever, we are cautious not to over-interpret the sign of the 
residual RM since the smooth, foreground RM surface 
subtraction will not have removed any small angular- 
scale variations in the RM of our Galaxy. Therefore, 
residual RM contamination from our Galaxy may con- 
tribute to local enhancements in the residual RM mag- 
nitude as well as changes in the sign of the residual RM 
in patches throughout the lobes. Therefore, we defer a 
detailed analysis of individual features of RM amplitude 
and sign to a more detailed study of the polarisation 
properties of the giant lobes. 

3.2. Origin of the Residual RM Signal 

The observed residual RM signal may be due to Fara- 
day rotation from thermal gas within the lobes, a thin 
boundary layer of swept up material surrounding the 
lobes or a large-scale fluctuation in the Galactic fore- 
ground Faraday depth coinciding with the position of 
the giant lobes. It is important to note that the subtrac- 
tion using the fit to the background RMs just outside 
the lobes not only removes the smooth Faraday depth 
contribution of our Galaxy but also of any smooth RM 
component on large angular scales caused by the intra- 
group medium or an extended halo of magnetized ther- 
mal plasma surrounding the lobes. 

The sign of the residual RM of the background sources 
within the lobes appears to correlate well with the RM 
from the diffuse emission (Figure [3]) . To quantitatively 
investigate the relation, we plot the residual background 
source RMs (RMbkg) against the corresponding resid- 
ual RMs of the polarized emission from Centaurus A 
(RMdiffuse) along the same sightline (Figure]!]). There are 
some obvious discrepancies, which are not unexpected 
given the different scales probed by the ATCA (~50") 
and Parkes (~16') observations as w ell as the intrinsi c 
scatter in the background source RMs (Feain et al. 2009 1 . 
A correlation-coefficient of 0.54 between these two quan- 
tities shows that both the sign and amplitude of the two 
different RM measurements are related, with a best fit 
to the data giving RM bkg ~ 1.3RM di g usc . 

Considering only the lines of sight in which the sign 
of RMbkg and RMdiffuse agree, we find a median ratio, 
q = RMbkg /RMdiffuse ~ 1-5. The median, rather than 
the mean, is quoted to avoid the results being skewed by 




Right Ascension (JSO0O) 

Fig. 3. — The residual RM (color scale ranging from —40 to 
+40 rad m~ 2 ) associated with Centaurus A after subtraction of the 
Faraday rotation due to foreground magnctoionic material. The 
residual RMs of the background sources inside the lobes are rep- 
resented by circles whose size corresponds to the magnitude of the 
RM while red/blue circles indicate positive/negative signed RMs. 
The legend in the top left corner gives the magnitude of the RM, 
in rad m — 2 , in relation to the size of the circles. The are no point 
source RMs towards the inner regions of the lobes because these 
regions were maske d out due to imag ing errors mainly caused by 
the bright core; see |Feain et all] ( |2009^ for details. 
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Fig. 4. — Plot of the background source residual RMs (RMbkg) 
versus the residual RMs from the polarized emission of Centaurus A 
(RMdiffuao) at the same position, with both axes in rad m — 2 . Solid 
line represents a best-fit line with slope 1.3 ± 0.1. 

a small number of outliers. Such a relation, with q > 1, 
is expected for most internal Faraday rotation models 
with q = 2 in the case of uniformly mixed e mitting and 
rotating regions (e.g., Cioffi & Jones 1980). However, 
a value of q ~ 2 is also expected in the case of a thin 
boundary layer of magnetoionic material surrounding the 
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Fig. 5. — Fractional polarisation image ranging from 0-40% 
overlaid by the total intensity contours. A particularly interesting 
feature is the high fractional polarisation located ~0.5° South- West 
of the core. The fractional polarisation is not shown out to the edge 
of the total intensity contours because a dramatic increase that is 
unrelated to the physical properties of Centaurus A is observed. 
The effect is as a result of the steep falloff in total intensity while 
the Galactic polarized emission remains approximately constant. 

lobes. This is assuming that the lobes and their imme- 
diate environments are symmetric in the sense that any 
postulated Faraday rotating region on the near side of the 
lobe is replicated on the far side. For the Northern lobe 
only, q ~ 1.3, while for the Southern lobe q ~ 1.5. This 
tentatively suggests that either a larger amount of gas 
is present within the Southern lobe or the lobe orienta- 
tion is such that the Southern lobe is further away; how- 
ever, it is difficult to claim this difference as significant 
given the smaller number of data points in the North- 
ern lobe compared to the Southern lobe. Overall, the 
range 1.3 < q < 1.5 supports an interpretation of in- 
ternal Faraday rotation in which there is an asymmetric 
distribution of emitting and rot ating regions within the 



lobes (e.g., Sokoloff et al. 1998), possibly with a small 



amount of Faraday rotation external to the lobes. How- 
ever, we cannot discount entirely external effects due to 
the orientation of the lobes or residual RM contamination 
from the Galaxy (as described in Section 3.1). We avoid 
any in-depth modelling based on the value of q due to the 
sparse and unevenly sampled nature of the background 
sources. A much denser grid of background sources is 
required for more robust comparisons and detailed mod- 
elling of the distribution of emitting and rotating regions 
associated with the giant lobes. 



skin surrounding the lobe boundary. If we assume that 
the residual large-scale RM signal we detect is due to a 
thin boundary layer of depth 6 ~ 20 kpc, then with an 
estimate of the line-of-sight magnetic field, Bu , we can 
obtain the expected electron number density of the skin, 
n e,skin- Based on the value of the magnetic field strength 
in the lobes of B ~ 0.9 ^G, as derived in |Abdo et al.| 



(2010) based on broadband modelling of radio and 7-ray 
data, we can place a strong upper limit for the line-of- 
sight magnetic field in the skin of B\\ ~ B/v3 ^0.5 /xG. 
This leads to a conservative lower limit of n e _ s ki n ~ 1.5 x 
10~ 3 (RM/12rad m~ 2 )(B eq /0.5 /iG) _1 (5/20 kpc 
This type of mod el is mainly based on the work of |Bick 



1 cm 3 . 



Previous work by Feain et al. ( 2009 1 found a turbulent 
RM signal associated with the southern lobe due to ei- 
ther turbulent structure throughout the lobe or in a thin 



nell et al.| ( 1990 1 who used hydrodynamic simulations 
to show that Kelvin-Helmholtz instabilities may form 
on the surface of radio lobes and potentially cause 
advection of the lobe magnetic field into the surrounding 
medium. 

The lobes of Centaurus A extend in to a sparse group 
of galaxies ( Karachentsev et al.||2007 ) with the Centau- 
rus A host galaxy, JNGC 5128 at the centre. An inter- 
galactic medium (I GM) density of r ti om ~ 10~ 3 cm~ 3 
was suggested by |Bouchard et al. (2007) to explain 
a discontinuity in the 111 properties ot Centaurus A 
group dwarf galaxies through ram-pressure stripping ar- 
guments. However, this would seem excessively large for 
such a poor group of galaxies. In order to estimate a more 
plausible intragroup density at the extent (r) of the giant 
lobes, we use a standard profile n igm (r) ~ n (r /a )~ b . 
This gives n iffm (200 kpc) ~ 10~ 4 cm~ 3 for the typi- 
cally derive d values of rip ~ 10~ 2 cm" 3 , q ~ 10 kpc 
and b ~ 1.5 ( |Mulchaey fc Zabludoff|[l998l ISunpOl^ , al- 
though the large scatter in the scaling parameters should 
be noted. The expansion of the lobes is expected to sweep 
up and possibly compress the intragroup medium around 
the edges of the lobes. Even if we consider a maximum 
compression factor of four, from a strong shock (which 
would be highly unlikely for the giant lobes), it is still not 
sufficient to reach n e ^ s ki n ~ 10 -3 cm -3 . Furthermore, 
a rather low magnetisation of the IGM is expected on 
scales of hundreds of kiloparsecs from the group centre, 
certainly much less than the equipartition value in the 
lobes, unless some additional processes are able to sub- 
stantially amplify th e IGM magnetic field, as discussed in 
Bicknell et al. ( 1990 ) for example. Many theoretical mod- 
els for the structure of extended lobes actually assume 
zero magnetisatio n of the ambient medium in to which the 
lobes evolve (e.g., Gourgouliatos et al.|2010[ ). Hence, the 
estimated value oin e , s kin ~ 10 -3 cm -3 should be consid- 
ered as a very conservative lower limit, corresponding to 
the maximum magnetisation of the IGM in the vicinity 
of the expanding lobes. Therefore, from the above argu- 
ments we conclude that the observed residual RM signal 
cannot be entirely explained by a thin skin scenario and 
that a significant fraction of thermal gas must be mixed 
throughout the lobe volume. 

3.3. Depolarization of the Lobes 

From the residual RM signal alone, we cannot deter- 
mine exactly how much thermal material is likely to be 
mixed in with the relativistic lobe plasma. A potential 
discriminant comes from analysing the change in frac- 
tional polarisation with wavelength. For mixed radio- 
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emitting and Faraday-rotating regions, the polarization 
angle of radiation produced at different depths within 
the lobes is rotated by different amounts which leads 



to wavelength-depe ndent depolarisation (Burn f966 
Sokoloff et al. f998). A thin-skin or boundary layer ef- 
fectively acts as a foreground Faraday sc reen and does 



not d epolarize the emission from the lobe ( Bicknell et al. 
19901. An unfortunate complication for the foreground 



scree ns is an effect known as external Faraday disper- 
sion (|Burn||1966 Tribble||1991|. This occurs when many 



turbulent ceils of magnetoionic material are within the 
telescope beam resulting in different amounts of Fara- 
day rotation along different lines of sight which causes 
depolarisation when averaged across the beam area. 

In order to better constrain the origin of the residual 
RM signal, we analysed the change in the degree of po- 
larisation between 2.3, 1.4 and 0.96 GHz. The 2.3 GHz 
image was obtained from the S-band Polarization All Sky 
Survey (S-PASS), which recovers the absolute level of 



the p olarized emission across the Southern sky ( Carretti 
2011 1 and hence provides a very robust estimate for the 



level of polarized emission from Centaurus A at 2.3 GHz. 
The 0.96 GHz fractional polarisation imag e was created 
using the Stokes / and p contour plots from | Cooper et ah] 
(1965). Digital copies of the se images ar e not a vailable 
so we used the algorithm of |Westphalen (1995) to reli- 
ably represent the contour plots in digital form. Given 
that the noise is not recovered in the digitized images, 
we limit our analysis to the highest signal-to - noise (S/N) 
areas where the estimated error is < 10% ( |Westphalen| 
1995 1. Both the 2.3 and 1.4 GHz images were smoothed, 
in Q and U to the resolution of the 0.96 GHz image. 

3.4. External Depolarisation Model 

For our depolarisation analysis we have chosen four 
positions which have high S/N as well as small observed 
variation in RM across the beam so that the depolarising 
effects of external Faraday dispersion are minimized. To 
model the depolarising effect of a foreground screen with 
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Fig. 6. — Plots of the degree of polarisation (p) versus wavelength 
squared (A 2 ) for four regions of the giant lobes of Centaurus A at 
2.3, 1.4 and 0.96 GHz. We chose areas with high S/N in p as 
well as small variations in RM across an area equivalent to the 
0.96 GHz beamsize. The dotted (red) curves represent a model 
of depolarisation solely from external Faraday dispersion while the 
dashed (blue) curves represent a model combining both internal 
Faraday rotation and external Faraday dispersion. In the top right 
and top left plots, the dashed (blue) and dotted (red) curves are 
almost identical. 



a Gaussian random magnetic field we use 
p = Po cxp(-2ctr M A 4 ), 



(3) 



where po is the intrinsic degree of polarized emission 
from the lobes and ctrm describes the RM fluctuations 
on scales smaller than the observing beam. 

Figure [6] shows the fractional polarisation versus wave- 
length squared for two positions in the Northern lobe, 
(13^27™, -42°29') and (13 h 26 m , -41°25'), and two 
in the Southern lobe, (13 ft 23 m , -43° 49') and (13 /l 24 m , 
— 44°21'). For the chosen positions, the variation in RM 
across an area the size of the 0.96 GHz beam is less than 
the estimated error of the individual RM measurements. 
The dotted (red) curves in Figure [6] represent a fit to the 
data using equation 3 with values of <7rm ranging from 
3.5-6.3 rad m~ 2 . Table 1 lists the values of the fitted 
parameters with their respective error estimates. From 
this we see that depolarisation from an external Faraday 
screen can, in principle, explain the total depolarisation 
observed at these four positions of the giant lobes. How- 
ever, we note that fluctuations in the complex polarisa- 
tion can become dominant at long wavelengths resulting 
in a much weaker than expected decrease in the d egree of 
polarisation due to extern al Faraday dispersion (Tribble 
1991| iSokoloff et al.||1998[ ). Hence, for Centaurus A, the 
observed depolarisation may not be very well described 
by an exponential decrease in the degree of polarisation. 

In Figure p] we can see a strong asymmetry between 
the polarizea intensity in the inner parts of the North- 
ern and Southern giant lobes, the origin of which could 



be attributed to the La ing-Garrington effect ( Laing 1988 
Garrington et al. 1988 ). This effect is commonly assignee 



to a halo of thermal gas, threaded by a t urbulent mag- 
netic field, sur rounding the radio galaxy (Garrington & 
Conway || 1 99 1[ ). In this case, the depolarisation is larger 
for the assumed more distant Southern lobe, in which the 
path length through the surrounding gas is larger. How- 
ever, there is no compelling evidence for the exact orien- 
tation of the lobes. The effect may also be explained by 
internal differences between the two lobes and/or from a 
gradient, with Galactic latitude, in the properties of the 
foreground Galactic Faraday screen. 

In order to comment further on this issue, we show 
the fractional polarisation within the giant lobes in Fig- 
ure 5| One of the most interesting features of this image 
is the high degree of polarisation located ^0.5° South- 
West of the bright core. We speculate that this region 
may be the oppositely-directed counterpart of the wel l 



studied Northern middle lobe (|Morganti et al 
which until now had no observe' 



1999 ) 

Southern counterpart 



This "Southern middle lobe" may have lost its structure 
and merged with the outer/older radio lobe after uplift- 
ing and mixing with a large amount of gas from the host 
galaxy. This scenario would require some asymmetry in 
the distribution of gas within the host galaxy whic h is not 
unrea sonable in the aftermath of a merger event ( Sparke 
19961 



IiTfact, a recent study by |Bell k, Comeau" (2012 1 found, 
for sources exhibiting the Lamg-Garrmgton effect, that 
the jet brightness asymmetry on kpc-scales cannot be 
explained by beaming and therefore must be intrinsic (in 
cases where the estimated kpc-scale outflow speeds are 
close to 0.1c). They conclude that a Laing-Garrington 
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TABLE 1 

External Faraday dispersion model 
parameters. 



Position 



(%) 



(rad m -2 ) 



(13 h 27 m , 


-42°29') 


39.3 ±0.4 


3.5 ±0.3 


(13 h 26 m , 


-41°25') 
-43°49') 


22.5 ±0.2 


5.5 ±0.3 


(13 A 23 m , 


34.0 ± 0.5 


4.2 ±0.4 


(13 h 24 m , 


-44°21') 


29.4 ± 0.6 


6.3 ±0.5 



Note. — Col. 1: Position of extracted data in 
RA and Dec in J2000. Col. 2: intrinsic degree of 
polarisation. Col. 3: variation in the external RM 
on scales smaller than the beam. 



TABLE 2 

Internal Faraday rotation model parameters. 



Position 


W 
(96) 


(rad m -2 ) 


CRM 
(rad m -2 ) 


(13 h 27 m , 
(13 ft 26 m , 
(13 h 23 m , 
(13 h 24 m , 


-42°29') 
-41°25') 
-43°49') 
-44°21') 


10.2 ± 1.7 

10.0 ± 0.7 

10.1 ± 1.2 
10.1 ± 0.9 


1.7 ± 0.7 

3.0 ±0.2 

2.1 ± 0.4 
2.1 ± 0.6 


1.2 ±0.6 
4.1 ± 0.3 
2.7 ± 0.4 
5.5 ±0.4 



Note. — Col. 1: Position of extracted data in RA and Dec 
in J2000. Col. 2: internal Faraday depth of the lobes. Col. 3: 
internal Faraday dispersion within the lobes. Col. 3: variation 
in the external RM on scales smaller than the beam. 

effect due to intrinsic differences in the Faraday rotating 
material within radio galaxy lobes cannot be discounted. 
Based on these results and our arguments in Section 3.2 
for a significant fraction of thermal material mixed in 
with the relativistic plasma of the lobes, we next consider 
a model in which a significant fraction of the observed 
depolarisation (and residual RM signal) is due to internal 
Faraday rotation. 

3.5. Internal Depolarisation Model 

If the observed Faraday rotation and depolarisation are 
due to thermal gas within the lobes, we can estimate 
the number density of the gas using an internal Faraday 
rotation model that accounts for depolarisation from the 
ordered magnetic field (B orc i) as well as from the random 
magnetic field (B r d m ) within the lobes. In this model, 
we consider the simplest case of a Gaussian distribution 
of Faraday depths with mean 4> and standard deviation 
as in a slab with a linear depth L along the line of sight 
(|Burn||1966[ |Sokoloff et aT|^998| . The intrinsic (A = 0) 
fraction of polarized emission from Centaurus A (po) is 
then modified following 



1 



-2crjA 4 +2iA^ 



P = Po 



-2ot 



2a 2 X 4 - 2i\ 2 <p 



(4) 



where o& — O.&lneBfd^Ld) 1 / 2 with d representing the 
scale of magnetic field fluctuations of the random mag- 
netic field (i?rdm) such that the number of cells along the 
line of sight is L/d. We have also included the effect of 
external Faraday dispersion, as described in Section 3.4, 
and use the same values of po . 

The dashed (blue) curves in Figure [6] represent fits to 
the data using the above model for each of the four lobe 
positions. In all cases, we find the magnitude of the in- 



ternal Faraday depth, \4>\ ~ 10 rad m 2 with o~& ranging 
from 1.7 to 3.0 rad m -2 and ctrm ranging from 1.2 to 
5.5 rad m -2 . Since we have only three data points and 
have to fit for three variables, this essentially guaran- 
tees a good fit, perhaps giving false confidence in the 
appropriateness of the model. Hence, we consider the 
fitted parameters of this model, shown in Table 2, as 
reasonable estimates of the model parameters which can 
describe the observed depolarisation given our current 
observational knowledge. 

Due to the random component of the magnetic field in 
the lobes, we expect to see an excess dispersion in the 
RM within the lobes over the RM dispersion outside the 
lobes (i.e., However, this is difficult to disentangle 
from the turbulent magnetoionic medium of our Galaxy 
that likely varies throughout the entire field as well as 
from effects more local to the source due to the orienta- 
tion of the lobes (e.g. the Laing-Garrington effect). The 
standard deviation of the RM from the diffuse Galactic 
emission (outside the lobes) is ~12 rad m -2 while the RM 
of the polarized emission from Centaurus A has a stan- 
dard deviation of ^14 rad m~ 2 . This difference cannot 
be regarded as significant, given our individual measure- 
ment errors (Section 2.3), but it is not inconsistent with 
the fitted values of a& used to describe the depolarisation 
in the lobes. 

The total magnetic field strength (B tot ) in the North- 
ern and Southern lobes has been estimated at ~ 0.9 /iG 



from the modelling of radio and 7-ray data ( Abdo et al. 



2010 1 . The intrinsic degree of polarization of synchrotron 
radiation (jpi) is reduced from its maximum level ac- 
cording to the relation p Q — PiB 2 rd j B 2 ot where B 2 ot = 

B lrd + B rdm an d Pi = ( 3 - 3a) /(5 - 3a), with a repre- 
senting the spectral index ( Burn|1966| . Therefore, using 
the spectral in dex of particular region s of the giant lobes 
obtained from Hardcastle et al. ( 2009 ) and the mean ob- 
served fractional polarization from our observations, we 
can estimate the strength of the ordered (B or d) and ran- 
dom (B rc im) magnetic field components. For example, 
at the position (13 h 23 m , -41°09') we have p = 22.5%, 
therefore B or d ~ B r d m ~ -By ~ 0.5 fiG. Considering 
a path length through the lobes of L ~ 200 kpc and 
a magnetic field fluctuation scale of d ~ 20 kpc, we 
use a& — (5 x 10 4 )n e i? I . dm and cj> ~ (2.5 x l0 5 )n e B\\ 
to find, for the four selected regions of the lobes, n e ~ 
( 0.9- 1.3) x 1Q- 4 cm - 3 . 
Feain et al. ( 2009[ ) provided an upper limit on the 



volume-averaged thermal electron density of the lobes 
of (n e ) < 1.3 x 10 _4 (0.5/B|i) _1 cm~ 3 , where in this case 
we use -B|| ~ 0.5 fiG instead of the equipartition value of 
1.3 fiG used in their paper. This limit was based on sam- 
pling of the background sources inside the lobes, probing 
a distribution presumably more inhomogeneous than the 
sparse sampling of the background sources could detect. 
However, the limit is not inconsistent with our results. 
Our measured value of n e ~ 10~ 4 cm -3 based on our 
internal Faraday rotation interpretation provides a very 
good description of all the current observables (i.e., the 
residual RM signal, the depolarisation and the excess 
RM dispersion associated with the lobes). Although we 
cannot definitively rule out external Faraday effects as 
an explanation for the total amount of observed depo- 
larisation, we prefer the internal Faraday rotation inter- 
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pretation since it better explains all the observational 
evidence. Taking this as our preferred model we now 
investigate some of the important consequences for the 
giant lobes. 

3.5.1. Implications for the Giant Lobes 

In order to calculate the pressure in the lobes due to the 
thermal gas (p t h ~ n e kT), we require an estimate for the 
temperature (T) of the gas. Recent X-ray results from a 
small region of the southern lobe found an excess of dif- 
fuse, thermal emission from Centaurus A over the back- 
ground, with a best-fit gas temperature o f kT ~ 0.5 keV 
and a number density of ~ 10~ 4 cm~ 3 ( jStawarz et al 



20121, in excellent agreement with the density estimate 
made in Section 3.3 above. Interestingly, this leads to 
a thermal pressure, p t h ~ 8 x 10~ 14 erg cm~ 3 which 
is approximately equal to the total non-t hermal pres- 
sure, (p e + + Ub) ~ 8 x 1 0~ 14 erg cm -3 (|Abdo et al. 



2010 Stawarz et al.| 2012|, where the magnetic energy 
density Ub = B^ ot /8n. It is important to note that 
our derivations do not depend on an equipartition as- 
sumption, and therefore, on the unknown contribution 
of relativistic protons. They depend instead on the mag- 
netic field value de rived from the broadband modelling in 



Abdo et al. ( 2010 ), whose primary assumption is that the 
detected 7-ray flux is solely due to the inverse-Compton 
emission of the radio emitting electrons. The derived 
magnetic field value is thus consistent with the case of an 
approximate energy equipartition between the magnetic 
field, radiating electrons, and potentially any relativistic 
protons present (i.e. Ub ~ U e > U p ). 

The majority of extended lobes in low-power radio 
galaxies are found to be under-pressured or in approx- 
imate pressure equilibrium with their external environ- 
ment, under the assumption of energy equipartition be- 
tween the magne tic field and radio-em i tting electrons 



withi n the lobes (Morganti et al. 1988 Croston et al. 
2008 1. Hence, those lobes must be either tar from the 
magnetic field-relativistic electrons energy equipartition, 
or there must be additional pressure provided by par- 
ticles contributing only weakly to the observed non- 
thermal continuum of the lobes to maintain the struc- 
ture of the expanding cavity. Here we have found direct 
evidence for thermal particles contributing substantially 
to the internal pressure in the giant lobes of Centau- 
rus A. Even so, approximate equipartition between the 
thermal and non-thermal gas pressures still cannot en- 
tirely explain the pressure mis-match in all low-power 
radio galaxies, indicating that relativistic protons may 
have a substantial contribution in some cases. 



Based on the derived number density of 
and the estimated volume V 



10- 



2 x 10 71 cm 3 



tie et al. 


2009 


Abdo et al. 


2010) 



(Hardcas- 



Hh 



nthmHfvV ~ 2 x 10 M , where m# is the mass of ion- 
ized hydrogen and we use a volume filling factor fy ~ 1 . 
Similar amounts of gas inferred from metal-enriched out- 
flows hav e also been observed in a number of other radio 
galaxies ([S imionescu et al. 2009 Kirkpatrick et al.|2009| 



McNamara & JNulscn 2012)! We consider the possibil- 
ity that the vast majority of the thermal material has 
either been entrained as the jet ploughed thr ough the 
interstellar m edium (ISM) of the host galaxy (Laing & 
|Bridle|[2002[ ) and/or pushed out by the over-pressured 



lobes as they expand outwards ( Begelman fc Cioffi|1989 
Churazov et al.pOOl ). Entrainment of thermal materia 
is considered as the most likely mec hanism for deceler 
ating relativistic jets (Bicknell 1994) but the amo unt of 
material requir ed is of the order of 10 -3 M Q /yr (Laing 
& Bridle 2002). Therefore, this is highly unlikely to be 



able to provide the necessary 200 M Q /yr requir ed by our 
results, for a lo be age of the order of 100 Myr (Hardcas- 
tle et al.][20 09j|. More likely is that the expanding lobes 
push out large amounts of gas from the host galaxy at- 
mosphere/halo through successive episodes of jet activ- 
ity. In Centaurus A there is evidence for at least four 
episodes of activ ity: the giant outer l obes, the North- 
ern middle lobe (|Morganti et al. 1999), the inn er lobes 
(IB urns et al. 1983p , and the parsec-scale jets (Tingay 
et al.|| 200l| ). 



Recent simulations (Wagner & Bicknell 2011) have 
shown the importance of inhomogeneity of the ISM which 
can cause deflection of the jet flow, allowing the radio 
source to effect a much larger volume of the host galaxy. 
In this way, a large amount of gas can be swept up by 
the influence of the jet at large distances from the center 
of the galaxy, possibly explaining the substantial amount 
of gas we have detected within the lobes. Such a large 
removal of gas from the host galaxy provides a direct 
mechanism for suppressing star formation while also lim- 
iting the amount of gas that can cool and fall back to- 
wards the center of the galaxy, which may also limit the 
growth of the central supermassive black hole. It is also 
possible that a significant fraction of the mass of thermal 
gas within the lobes was entrained from the intergroup 
medium and dispersed over the entire lobe volume. More 
detailed simulations with synthetic Faraday rotation ob- 
servations will be required to help distinguish between 
these scenarios. 

4. CONCLUSIONS 

We have presented results from a spectropolarimetric 
study, using the Parkes radio telescope at 1.4 GHz, in- 
vestigating the Faraday rotation of the diffuse polarized 
emission from the giant lobes of the radio galaxy, Cen- 
taurus A. Using previous results from an RM g rid of 
background radio sources just outside the lobes (Fcain 
et al. 2009), the smooth foreground contribution to the 
observed RM from Centaurus A was subtracted, leaving 
a mean residual RM signal of ~ 12 rad m~ 2 . 

Investigation of whether the residual RM signal comes 
from a thin-skin/boundary layer of magnetoionic mate- 
rial surrounding the lobes or from thermal gas internal 
to the lobes, found the thin-skin scenario highly unlikely 
given that it requires at least an order of magnitude en- 
hancement of the swept up gas over the expected intra- 
group density on these scales. It should be noted that 
we cannot conclusively rule out that the residual RM sig- 
nal comes from a degree-scale fluctuation in the smooth 
Galactic foreground that was not sampled by the ensem- 
ble of background sources outside the lobes and happens 
to align with the orientation of the giant lobes. 

From our investigation of the degree of polarisation 
at three separate frequencies (2.3, 1.4 & 0.96 GHz) at 
four positions of high S/N, we find strong depolarisation 
which cannot be explained solely by the effects of beam 
depolarisation. Considering all the available data, we 
consider depolarisation due to Faraday rotation within 
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the lobes as our preferred scenario, even in the presence 
of external Faraday dispersion from a foreground screen. 
In this case, we estimate the number density of the ther- 
mal gas n e ~ 1CP 4 cm~ 3 which gives a total ionized gas 
mass within the lobes of ~ 10 10 M Q . Recent X-ray ob- 
servations of a small region of the southern lobe also find 
evidence for therma l gas within the lobe with a temper- 
ature of ~ 0.5 keV (Stawarz et al. 2012). From this we 
estimate that the thermal pressure within the lobes is 
approximately equal to the non-thermal pressure imply- 
ing that the thermal gas, radio-emitting electrons and 
magnetic held are all in approximate equipartition with 
each other. 

Future radio polarisation observations are required to 
have much wider A 2 -coverage to provide the high resolu- 
tion in Faraday depth space needed to uniquely separate 
the internal and external Faraday rotating components. 
Indeed the planned Global Magneto-Ionic Me dium Sur 



vey (GMIM S) which covers 300-1800 MHz ( |Wolleben 
et al.|[2009 1 provides an excellent opportunity to resolve 



the Faraday rotation structure in the giant lobes of Cen- 
taurus A. 
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APPENDIX 



For observations which have more than one source of polarized emission along the line of sight with different Faraday 
depths and/or Faraday thick nesses, the measured RM can vary with wavelength (Law et al. 2011 Farnsworth et al 



2011 O'Sullivan et al. 2012). Since we detect both the diffuse polarized emission from our Galaxy as well as the 



polarized emission from Centaurus A, we need to consider how this may affect our measured RM. The RM resolution 
of 310 rad m~ 2 limits our ability to distinguish between multiple regions of polarized emission whose Faraday depth 
differs by less than this amount. Therefore, we use models of two line-of-sight polarized components with different 
Faraday depths to investigate at what points within the lobes of Centaurus A our measured RMs accurately represent 
the true Faraday depth of the polarized emission from Centaurus A. 
To illustrate the effect, we define a two component model 



1 



Po 



2al A 4 - 2i\ 2 c 



,2i( X o+<t>e* 



- Vg 



<M 2 



(5) 



in which the polarized emission of our Galaxy (p g ) is described by a uniform slab (Burn 1966) of Faraday depth <j) g 
while the Centaurus A model is described in Section 3.5. We use the values for position (13' i 27 m , — 42°29') listed in 
Table 2 with the Galactic Faraday depth, tp g = —58 rad m~ 2 taken from the fit for the foreground RM surface. The 
intrinsic polarization angle of the Centaurus A emission, xo — 30° , is taken from a derotated A = polarization angle 
at this position. For the Galactic emission, we use a completely arbitrary value of xog = 90°. We then set po = 1.0 
and plot the results for various fractions oip g /po, shown in Figure 7. 

The top panel of Figure 7 shows the Faraday rotation measure calculated as RM = d\j dX 2 , where X—\ arctan U /Q 
and P = Q + ill. This shows that as long as p g /po < 0.15 the possible systematic error in the measured RM is less 
than 3 rad m -2 , where the expected RM is \<j) + 4>g (solid horizontal line). As the Galactic polarized emission (mean 
value ~52 mJy) becomes a larger fraction of the Centaurus A polarized emission, the systematic error in the measured 
RM becomes much larger (dotted line). Since we are unable to disentangle the two contributions with our current 
dataset, we have only included in our calculations the regions within the lobes of Centaurus A where the polarized 
emission is > 300 mJy beam . This value roughly corresponds to p g /po < 0.15, which means that even in the regions 
of lowest polarized intensity the potential systematic error in the RM is less the mean measurement error within the 
lobes of ~3 rad m~ 2 . 
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Fig. 7. — Plot of RM vs. A 2 , to illustrate the RM variation due to sampling multiple regions of mixed polarized emission and Faraday 
rotation along a single line of sight. The RM is calculated as dx/d\ 2 and the solid horizontal line is the expected RM (see Appendix for 



details). Solid line: p g /p c = 
covered by our observations. 



0.15, dashed line: p g /p c = 0.3, dotted line: p g /p c = 0.5. The shaded areas indicate the range of A space 



